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ABSTRACT 

We  consider  a  one-dimensional  model  problem  for  the  motion  of  a  viscoelastic  material 
with  fading  memory  governed  by  a  quasilinear  hyperbolic  system  of  integrodifferential 
equations  of  Volterra  type.  For  given  Cauchy  data  in  Z,°°(R)  n  L2(R),  we  use  the  method 
of  vanishing  viscosity  and  techniques  of  compensated  compactness  to  obtain  the  existence 
of  a  weak  solution  (in  the  class  of  bounded  measurable  functions)  in  a  special  case. 
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1.  INTRODUCTION. 


COPY  I 
lNSPtCT£oj 


The  purpose  of  this 

paper  is 

global  weak  solutions,  in 

the  class 

lem 

( 

wt  =  vx, 

(VE) 

Vf  =  <TX, 

i 

tu(x,0)  = 

xeR,t>o, 

( VE )  <  =  **» 

{  te(i, 0)  =  w0(x),  v(x,0)  =  u0(*),  x  €  R, 

where  the  function  cr(x,t)  is  determined  by  the  history  of  u?(x,-)  through  the  constitutive 
assumption 


(CA) 


r(x,t)  =  if(w(x,t))  +  [  k(t  —  T )ifr(w(x,  r))di 

Jo 


The  given  functions  V’(w)  and  k(t)  are  assumed  to  be  smooth  and,  in  addition, 

>0,  w  G  R, 

so  that  the  structure  of  (VE)  is  hyperbolic.  In  this  paper  the  aforementioned  goal  is 
achieved  in  the  important  special  case  ib  =  v7?  when  the  data  tto»vo  €  Z/°°(R)  n  L2(R), 
using  the  method  of  vanishing  viscosity  and  techniques  of  compensated  compactness  (cf. 
[21,22,23,7,17]). 
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The  system  (YE)  is  a  model  for  one-dimensional  motion  of  a  viscoelastic  material  with 
fading  memory.  The  functions  v(x,t),iv(x,t)  and  cr(z,t)  stand  for  the  velocity,  deforma¬ 
tion  gradient  and  stress,  respectively,  while  the  constitutive  assumption  (CA)  states  that 
the  stress  is  a  particular  functional  of  the  history  of  the  deformation  gradient;  in  (CA)  the 
history  of  is  assumed  to  be  zero  for  t  <  0  and  the  body  force  is  taken  to  be  zero. 

Under  appropriate  assumptions  on  the  kernel  k,  the  memory  induces  a  weak  dissipative 
mechanism  whicn  competes  with  the  hyperbolic  character  of  (YE).  Under  such  assump¬ 
tions,  qualitative  properties  that  would  make  (YE)  a  reasonable  model  for  the  motion  of 
viscoelastic  materials  include: 

(a)  For  smooth  and  small  data  w0,vq ,  (VE)  should  possess  globally  defined  classical  so¬ 
lutions,  which  decay  to  equilibrium  as  t  —>  oo.  Such  results  have  been  established  in 
[12,5,20:  if  v  ~  p,  and  in  [6,9]  if  ip  ^  p;  for  more  general  models,  see  [18] . 

(b)  By  contrast,  if  the  smooth  data  are  chosen  sufficiently  large,  smooth  solutions  of  (VE) 
should  develop  singularities  in  their  derivatives  in  finite  time.  This  has  been  shown  in 
[2,18j;  for  other  models  see  [18]. 

(c)  For  arbitrary  L°°  data  wq,vq,  weak  solutions  should  exist  in  Rx  [0;oo).  This  property 
is  established  here  in  the  special  case  ip  =  p. 

There  is  a  parallel  theory  for  the  single  conservation  law  with  memory 


(CLM) 


wt  +  orx  =  0, 
u?(x,0)  =  u'o(x) 


where  a(x,t )  is  given  by  (CA).  For  this  equation  property  (c)  was  recently  established  by 
Dafer mos  [4;  using  the  method  of  compensated  compactness;  unfortunately,  this  elegant 
approach  does  not  seem  to  extend  to  (VE). 

In  the  sequel,  we  restrict  attention  to  the  case  p  =  iJj.  Then  (VE)  reads 


(1.1) 


wt  =  vx, 

vt  =  p(w)x  +  f*  k(t  ~  r)p(w(x,  r))xdr 
ic(x,0)  =  ie0(z),  v(x,0)  =  u0(x), 


i£R,  t  >  0, 
if  R. 


Regarding  properties  (a),  (b)  and  (c),  it  is  instructive  to  compare  (1.1)  with  the  equations 
of  one-dimensional  elasticity 


(£) 


rwt  =  vx 
l  Vt  =  <f(w)x, 


and  with  the  equations  of  one-dimensional  frictionally  damped  elastic  materials 


(FE) 


r  w(  -  vx 
\  Vt  = 


w)z  -  V. 


The  equations  of  elasticity  enjoy  property  (c)  (DiPerna  7 ] );  by  contrast,  property  (a)  is 
not  pertinent  to  (E),  since  smooth  solutions  of  (E)  generally  break  down  in  finite  time.  On 


2 


the  other  hand  (FE)  exhibits  a  closer  resemblance  to  (1.1)  enjoying  all  the  properties  (a), 
(b)  and  (c)  (see  ;]5,19,  7’  respectively). 

The  similarity  between  (1.1)  and  (FE)  is  revealed  by  the  following  transformation  due 
to  McCamv  [12;.  Let  r(t)  be  the  resolvent  kernel  associated  with  fc;  i.e.,  r  is  the  solution 
of  the  linear  Volterra  equation 

(r)  r(t)  +  f  k{t  —  T)r(T)dr  =  k(t),  t  >  0. 

Jo 

Convolving  (1.1)2  with  r(t),  a  simple  calculation  yields 


(1.2) 


/  k(t  -  t)<p(w(x,t))x<It  =  /  r(t  -  r)ut(i,r)dr 

Jo  Jo 

=  r(0)u(x,t)  —  r(t)v0(x)  -f  f  r'(t  —  r)u(x,r)d7 

Jo 


Thus,  for  smooth  solutions  (1.1)  is  equivalent  to 

ifR,  t  >  0, 


W,  =  V, 


(1.3) 


<  vt  =  <p(w)x  +  J’H, 
w(x,0)  =  w o(z),  v(x,0)  =  v0(x )  x  G  R, 


where 

(1.4) 


[v](r,  t)  :=  r(0)u(x,  t)  ~  r(t)vo(x)  +  f  r'(t  —  r)r(x,r)dr. 

Jo 


Under  physically  reasonable  assumptions  regarding  the  kernel  k  (e.g.  k  =  a',  with  a 
smooth,  positive,  decreasing  and  convex  on  [0, 00)),  r(0)  =  a'(0)  <  0  and  the  term 
r(0)v(x,t )  has  a  damping  effect.  The  effect  of  this  term  is  dominant  close  to  equilib¬ 
rium  thereby  inducing  property  (a).  However,  far  from  equilibrium,  the  hyperbolic  part 
of  (1.1)  dominates,  irrespective  of  the  sign  properties  of  k.  It  is  thus  conceivable  that  the 
analysis  of  weak  solutions  is  not  affected  by  the  sign  properties  of  the  kernel  k(t). 

We  assume:  the  constitutive  function  ip  satisfies 


(1.5) 


tp  :  R  — >  R  is  a  twice  continuously  differentiable 
function  such  that  <p'(w )  >  0,  w  6  R; 
p  has  a  single  inflection  point  at  w  =  uq  and  is 
convex  on  (uq,oo)  and  concave  on  (-00, uq). 


The  kernel  k  satisfies 


(1.6) 


k  :  [0, oc)  — ♦  R,  fcGC1[0,oc), 


and  the  data  il'o(x),  i’o(j’)  satisfy 


(1.7)  w0(x),  v0(x)  e  L°°( R)  n  L2( R). 

We  prove: 

Theorem  1.1.  Let  the  hypotheses  (1.5)-(1. 7)  be  satisfied.  Given  T  >  0,  there  exists  a 
weak  solution  {u(x ,t),v(x,t)}  of  ( 1.1)  on  R  x  [0,  T),  such  that 

(w,v)  e  L°°([0,  T];  L2(R))  n  L~(R  x  [0,T]). 

The  proof  is  carried  out  in  Sections  2,  3  and  4  using  the  method  of  compensated 
compactness  of  Murat  [13]  and  Tartar  [21,22,231.  This  approach  has  been  employed  with 
success  by  Tartar  i22]  to  obtain  Z,°°  solutions  for  the  general,  scalar  Burgers  equation,  by 
DiPcrna  17]  and  Rascle  :17  ,  to  construct  L°°  solutions  to  the  hyperbolic  system  (E),  and 
by  Dafermos  i4',  to  obtain  L°°  solutions  for  (CLM). 

The  weak  solutions  of  (1.1)  will  be  sought  as  the  e  \  0  limits  of  solutions  of  the 
regularized  problem 

'  wt  =  vx  +  ewxx, 

Vt  =  <p(w)z  +  Jo  k(t  -  r)ip(w(x,r))xdT  ( x,t )  6  R  x  (0 ,T], 

C1-8)  '  +e(u**  +  Jo  k{t  -  r)vIX(x,r)dr), 

.  iy(a:,0)  =  iu0(i),  u(x,0)  =  v0[x),  iER. 

This  unconventional  regularization  preserves  one  of  the  main  features  of  (1.1);  namely,  that 
the  integration  with  respect  to  t  offsets  differentiation  with  respect  to  x  and  the  memory 
terms  are  of  lower  order.  In  Section  2  the  system  (1.8)  is  transformed  into  a  parabolic 
regularization  of  (1.3).  It  is  also  proved  there  that  the  initial  value  problems  (1.1)  and 
(1.3)  are  equivalent  for  weak  solutions. 

In  Section  3  we  prove  an  existence  theorem  for  (1.8)  and  obtain  the  necessary  a-priori 
estimates  in  order  to  pass  to  the  limit  as  e  \  0.  Most  notably,  using  a  class  of  exponentially 
growing  convex  entropies  of  (E)  constructed  by  Dafermos  [3],  we  show’  that  the  solutions 
{wct'e}  of  (1.8)  lie  in  a  bounded  set  in  L°° ,  uniformly  in  e.  (In  view  of  the  nonlocal  nature 
of  the  memory  term,  it  does  not  appear  possible  to  establish  Z/°°-estimates  for  solutions  of 

(1.8)  by  constructing  invariant  regions.) 

The  results  of  Section  3  imply  that  a  subsequence  can  be  extracted  converging  in  L°° 
weak  star.  Under  the  L°°  weak  star  convergence,  composite  weak  limits  are  characterized 
as  expected  values  of  a  family  of  probability  measures  i/(x,f),  called  Young  measures 
(Tartar  '21,22,23  ).  In  Section  4  we  employ  the  techniques  of  DiPerna  [7]  concerning  the 
convergence  of  approximate  solutions  to  the  equations  of  elasticity  (E)  and  a  lemma  of 
Murat  14]  to  show  that  the  Young  measure  i^(x,<)  reduces  to  a  Dirac  mass  for  almost 
all  (x,f).  Thus,  a  subsequence  {u-V,vei}  of  solutions  of  (1.8)  converges  almost  everywhere 
and  the  limit  is  a  weak  solution  of  (1.1). 


Certain  results  concerning  weak  solutions  of  (VE)  have  appeared  in  the  literature. 
Greenberg  [8]  (also  see  18],  Section  2.6)  established  the  existence  of  travelling  wave  so¬ 
lutions  (steady  compression  shocks)  for  a  history  valued  problem  associated  with  (VF.). 
Boldrini  [lj,  used  the  compensated  compactness  method,  to  show  that  as  the  memory 
weakens  ( k  —  k{6,t)  =  Ci(£)),  uniformly  bounded  solutions  of  (1.8)  converge,  as 

6  1  0,  to  a  solution  of  the  elastic  problem  (E). 

2.  WEAK  SOLUTIONS;  THE  REGULARIZED  PROBLEM. 

We  seek  weak  solutions  of  (1.1)  in  the  space  L°°  of  bounded  measurable  functions. 
Within  this  framework,  a  pair  of  functions  {u)(x,  <),  v(x,  /)}  in  L°°(;0,  T)\  L2(R)  P  Z,°°(R)) 
is  a  weak  solution  of  (1.1)  on  [0 ,T]  x  R  if  it  satisfies 


(2.1) 


r 


oo 


oo 


£(x,  0)u>o(x)d;r  =  0, 


(2.2) 


/  /  [O  -  (x(v>(*>)  +  f  k[t  -  T)<p(w(-,T))dT)]dxdt 

J  0  J  oo  J  0 

+  f  ((x,0)v0(x)dx  =  0, 

J  —  OO 


for  any  pair  {£(z,  <),  £(x.  <)}  of  C1  functions  with  compact  support  on  R  x  [0,  T]  satisfying 

«*,r)  =  C(x,r)  =  o. 

The  weak  solutions  of  (1.1)  will  be  constructed  as  the  e  \  0  limits  of  solutions  of 
the  regularized  system  (1.8).  The  relevance  of  the  unconventional  regularization  (1.8)  is 
revealed  by  the  following  calculation.  Convolving  (1.8)2  with  r(t),  the  resolvent  kernel  of 
k(t),  and  using  (r)  and  (1.4)  we  arrive  at 

/  k(t  -  t)(<p(w(x,  t))x  +  evxx(x,  r))dxd.T  —  /  r(t  -  r)vt[x,  r)dr 
Jo  Jo 

(2.3)  =  r(0)v(x,t)  -  r(t)v0{x)  +  f  r'(t  -  t)v(x,  r)dr  =  T[v](x,  t). 

Jo 

Thus,  the  initial  value  problem  (1.8)  can  be  written  in  the  form 


(2.4) 


(x,<)  6  R  x  (0,T] 


Wt  =  Vx  +  £WXX, 

vt  =  v{w)  x  +  T\y)  +  evxx, 
w(x,0)  —  u’o(x),  v(x,0)  =  n0(x),  x  6  R. 


In  performing  the  above  transformation,  a  technical  difficulty  arises  associated  with 
the  fact  that  the  initial  data  w0,  r0  are  only  assumed  to  be  in  L°°  P  L2 .  Experience  with 
parabolic  equations  suggests  that  we  should  not,  in  general,  expect  r„(x,-)  and  r((x,  ) 


to  be  integrable  on  [0,<]  for  solutions  of  (2.4).  Thus  we  need  to  make  (2.3)  precise.  It 
turns  out  that  (2.3)  holds  in  a  weak  sense.  Moreover,  (1.8)  is  equivalent  to  (2.4),  and 
furthermore,  (1.1)  is  equivalent  to  (1.3)  for  weak  solutions.  Specifically,  we  prove 

Proposition  2.1.  Let  u’o(.r),  v0(x)  6  L°°(R)  n  L2(R),  k(t)  £  C1 1 0,oo). 

(a)  The  initial  value  problems  (1.1)  and  (1.3)  are  equivalent  for  weak  solutions  w,v  of 
class  i“(;0,ri;/=°(R)nf2(R)). 

(b)  The  initial  value  problems  (1-8)  and  (2.4)  are  equivalent  for  weak  solutions  w,v  of  class 

L°°([ 0,T];  L°°(R)  -  L2( R))  n  l2([0,r;; 

PROOF.  We  first  consider  part  (a).  Let  {u>(x,  t),  v(x,  t )}  be  bounded  measurable  functions 

satisfying  (2.1)  and  (2.2).  We  proceed  to  show  that 

{u>(a\  t),  u(x,t)}isa  weak  solution  of  (1.3);  i.e.,  it  satisfies  (2.1)  and 

(2.5)  [  f  [T)tv  -  vMw)  +  r]T\v\)dxdt  +  /  q(x,  0)v0(x)dx  =  0 

Jo  J  oe  J-ao 

for  any  pair  {£(x,  t),  rj(x,  *)}  of  C1  functions  with  compact  support  on 
Rx  [0,  T],  satisfying  £(x,T)  =  r,(x,T )  =  0. 

Indeed,  given  any  test  function  Tj(xyt)  we  define 

(2.6)  C(x,  t)  =  v{x,t)  -  J  r{r  -  t)r1{x,T)dT. 

Then  ((x,t)  is  a  C1  function  with  compact  support  on  Rx[0,T]  and 

((x,T)  =  0.  Thus,  C(x,t)  can  be  used  as  a  test  function  in  (2.2).  With  £(*,/)  defined  by 

(2.6) ,  an  easy  computation  using  (r)  yields 

/  /  Cz{<p[w)+  f  k(t  -  T)ip(w(-,T))dT)dxdt 

J  0  J  —  oo  Jo 

(2.7)  =  f  (  rjx‘fi(w)dxdt, 

J  0  J  —  oo 


and,  by  (1-4), 


noo  i*00 

C,tv  dxdt  +  /  ((x,  0)v0(x)dx  = 

•  oo  *  —  oo 

=f/:  ( Vt  +  r(0)rj  +  /  r'(r  -  t)q(-,T)dT)v  dxdt 

r(r)ri(x,T)dT^vo(x)dx 

pT  /•«= 

-  /  /  (r}tv  -f  r]T{v\f)dxdt  +  /  77(i,0)vo(*)dx. 


Combining  (2.2)  with  (2.7)  and  (2.8)  we  conclude  that  {u<(x,t),  v(x,  t)}  satisfies  (2.5)  and 
is  thus  a  weak  solution  of  (1.3). 

Conversely,  suppose  that  {w(x,  t).  r(x,  t)}  satisfies  (2.1)  and  (2.5).  Given  (( x,t ),  we 
solve  the  integral  equation  (2.6)  for  T)(x,t )  and  obtain 


T]{x,t)  =  ((x,t)  +  J  k(r  —  t)C(x,  r)dr. 


Then,  we  employ  q(x,t)  as  a  test  function  in  (2.5)  and  use  (2.7)  and  (2.8)  to  arrive  at 
(2.2).  This  concludes  the  proof  of  part  (a). 

By  comparing  the  relations  (2.2)  and  (2.5)  we  deduce  that 

(2.10)  f  f  \(F[v]  +  Cxf  k(t  -  t)v(w(-,t))(1t  dxdt  =  0, 

JO  J  —  oo  ^  Jo 

for  every  test  function  £(x,t).  This  is  a  weak  form  of  relation  (1.2). 

The  proof  of  part  (b)  follows  the  same  pattern  and  will  be  omitted.  Moreover,  the 
analogue  of  (2.10)  is  the  weak  form  of  (2.3). 

3.  A  PRIORI  ESTIMATES. 

Our  object  is  to  obtain  a-priori  estimates,  independent  of  e,  for  solutions  {u/,!;1}  of 
(1.8);  we  employ  the  equivalent  system  (2.4).  For  the  remainder  of  this  section  we  drop 
the  superscript  e. 

Let  {w(x,t).v(x,t)}  be  a  solution  of  (2.4),  corresponding  to  initial  data  iu0,t’o  in 
L°°nL2,  and  tending  to  (0,0)  as  |i|  — »  oo.  The  solution  is  assumed  sufficiently  smooth  to 
justify  the  steps  in  the  derivation  of  the  a-priori  estimates  (e.g.  the  regularity  properties  in 
Theorem  3.1  below  suffice).  In  the  sequel  C  will  stand  for  a  generic  constant  depending  on 
the  X°°(R)  and  L2(R)-norms  of  the  initial  data,  on  the  C1  [0,  Tl-norm  of  r(<),  on  properties 
of  the  function  <p(u>),  on  T  but  not  on  e.  Whenever  the  constant  depends  on  e  it  will  be 
denoted  by  Ce. 

Our  analysis  of  (2.4)  will  be  based  on  the  concept  of  entropy-entropy  flux  pairs  for  the 
elastic  problem  (E)  (cf.  Lax  [ll]).  A  smooth,  convex  function  tj(w,v)  defined  on  RxR  is 
an  entropy  for  (E),  with  corresponding  entropy  flux  q(w,v),  if 

(3.1)  dtT]{w{x,t),v(x}t))  4-  dxq{w(x,t),v(x,t))  =  0 

for  any  smooth  solution  {u’(  j,  t),  v(x,  f)}  of  (E).  Such  pairs  are  generated  as  solutions  of 
the  system  of  equations 

/o  ^  J  =  -v'WVv 


provided  rj(w,v)  is  convex.  Eliminating  q(w,v)  in  (3.2),  we  find  that  tj(w,v)  must  be  a 
convex  solution  of  the  linear  wave  equation 

(3.3)  Vww  = 


OT.TWii 


-Kwcy- 


q(w,  r)  is  then  determined  by  (3.2).  A  classical  example  of  an  entropy-entropy  flux  pair  is 


(3.4) 


V( 


1 

w,v)  =  -v2  +  J  Y>(()d(,  q(w,v)  =  -vifi(w). 


For  each  entropy-entropy  flux  pair  {r)(w,  v),  q(w,  r)}  the  functions 
{fj(w,v),q(iv,v)s  defined  by 


(3.5) 


f  fj(w,v)  :=  77(11;, r)  -  77(0,0)  -  77^(0,0)14'  -  q„(0,0)v, 

1  g(^,u)  :=  <7(777,7;)  -  9(0,0)  -f  77^(0,  0)u  -  T7„(0,  0)(^j(u’)  -  s?(0)), 


form  an  entropy-entropy  flux  pair  for  (E)  that  vanishes  to  quadratic  order  in  (tjlw)  at 
(0,0).  Along  solutions  {tt;(x,  t),  v(x,  t)}  of  (2.4),  a  simple  computation  yields  the  identity 


(3.6) 


dtr)(w,v)  +  dzq(  w,v)  =  fjv(w  ,v)T[v]  +  ed2fj(w,v) 

-  e[rjww(w,  v) w'l  +  2tju.v(w,  v)wxv.^  ~  Vw{w,v)i )\}. 


Integrating  (3.6)  over  Rx(0,<),  0  <  t  <  T,  we  obtain 

,00 

I  77(777(2;,  t),  v(x,  t))dx 

J  —OO 

+  e  /  {Vu=w{w,v)w2x  +  2t)wv(w,v)wxvx  + -qvv(w,v)vl) 

Jo  J~  oo 

/oo  ft  foo 

9(7Z70  (x),  770  (x))dx -f  /  /  f)v{w,v)T[v\dxdr,  0  <t<T. 

-oo  Jo  J  —  oo 


For  the  first  estimate,  we  employ  the  pair  (3.4)  and  use  (3.5)  and  (3.7)  to  deduce 

ft  *00 


[  [\v2(x,t)  +  $(w(x,i))}dx  +  e  f  f  (v'(w)wl  +  v2x)dxdr 

2-oo  4  Jo  J  —  00 

=  J  [\vo(x)  +  ^(w0(*))]^a! 

+  f  f  v!F[v}dxdT,  0  <  t  <  T, 

Jo  2-00 


(3.8) 
where 

(3.9)  $(w):=  f  (<f(0  -  ¥>{0))d£. 

Jo 

On  account  of  (1.5), 


(3.10) 


<p'(w)  >  >  0,  for  every  rv  6  R 


(3.11) 


$(w)  >  — — —  U’2  >  0,  for  every  w  £  R. 


Using  (1.4)  and  the  Cauchy-Schwarz  inequality, 


t  1*00 


a 


vTv  dxdr 


=  /  t>(i,r)  r(0)t>(x,  r)  -  r(r)r0(x)  4  /  r'(r  -  s)r(x,  s)ds\ dxdr 

J  0  J  —  oo  Jo  J 

a:  v2  dxdr^j 


(3.12)  <(?(l4-  //.:  r2(x,  r)dxdr j ,  0  <  t  < 

Combining  (3.8)  with  (3.10)  and  (3.12) 


l  . 


(3.13) 


/OO  A  OO 

v2(x,t)dx  <0  +  0  I  v2dxdr , 

■  oo  J  0  J  —  oo 


whence  by  the  Gronwall  inequality,  (3.8),  (3.10)  and  (3.11), 


/: 


[r2(x,*)  -f  w2(x,t)]dx 

ft  fOO 


(3.14) 


+  e 


j  j < 


w\  4-  v2x)dxdr  <  C,  0  <  t  <  T, 


where  C  is  independent  of  e. 

Our  next  objective  is  to  establish  a-priori  L°°  estimates,  independent  of  e,  for  solutions 
{u’(x,  <),  v(x,  <)}  of  (2.4).  We  follow  the  development  of  Dafermos  [3] .  The  following  facts 
are  proved  in  ; 3,  Section  2j:  For  t?(w)  as  in  (1.5),  the  wave  equation  (3.3)  admits  a  class 
of  solutions  {r/fc)(u',  r)}fc>0  on  RxR  which  are  strictly  convex  and  grow  exponentially  at 
infinity.  These  solutions  have  the  form 

(3.15)  t) ^k\w,v)  =  y^(u;)cosh  Ur,  0  <  k  <  oo, 
where  Y^k\w)  is  the  solution  of  the  initial  value  problem 

Y^y'{xv)  =  UV(^)V(fc)(u-), 

(3.16) 

r(*>(tr.)  =  1,  r<*>'(uq)  =  0,  o  <  Jfc  <  OO. 


The  functions  Y'k'(iv)  satisfy  the  estimates 

(3.17)  1  ^(ir)  ~  cosh  k y/p  '(u’,)(u;  -  u^)j,  — oo  <  w  <  oo,  0  <  k  <  oo, 

(3.18)  \Y<k)  (u;)j  <  ky^p\w)Y^k\w),  -oo  <  w  <  oo,  0  <  k  <  oc 
and 

»ll> 

(3.19)  Yr)(u’)<  exp;  A:  /  ■)/'•?' (J)d£] ,  -oo  <  w  <  oo,  0  <  A:  <  oc. 

J  tit, 

We  will  estimate  the  solution  {u.’(x,  t ),  r(x,  <)}  of  (2.4)  by  monitoring  the  evolution  of 

f  ij{k)(w{x,t),v(x,t))dx. 

J  —  OO 

In  view  of  the  convexity  of  fj(k\w,v),  (3.7)  yields 

/OO  /•  OO 

fj{k\w(x,t),v(x,t))dx  <  /  i7(lc)(u-0(i),uo(x))<ix 

oo  J  —  oo 

(3.20)  +  f  f  fj[k\w(x,T),v(x,T))!Flv)(x,T)dxdT,  0  <  i  <  T, 

J  0  J  —  oo 

where,  by  (3.5),  (3.15), 

(3.21)  =  77(fc)(™>”)  -  K(t)(0)  -  Y{k)'(0)w  >  0, 


(3.22)  fj[k\w,v)  =  A;(tanh  kv)^k\w,v). 

Using  (3.22),  (3.21),  (3.18),  the  Cauchy-Schwarz  inequality  and  (1.4)  the  last  integral  in 
(3.20)  can  be  estimated  as  follows: 


f}[k\w,  v)T[v)dxdT 


(3.23) 


=  /  J  A:(tanh  kv)  ^fj^k\w,v)  4-  Y^(0) 

+  {0)w^!F[v}dxdT 
<  Ck2Y^k\$)  f  f  ( [it? [  +  juDl^fv]! dxdr 

Jo  J -oo 

+  k  f  f  fj^-k\w,v)\J:'[v}]dxdT 

Jo  J  —  oo 

<Ck2Y^k\0)[f  f  (w2  +  v2)dxdr  -r  f  t'Q(x)dx] 

Jo  J -oo  J  —  oo 

+  Ck  j  M(t)[[  f)^k\w(x,  r),  v(x,  r))dx]dr, 

Jo  J- oo 


(3.24)  A/(<)  :=  sup{l  -fi  jv(i,  t)\  +  /  |t>(x,  r)i<fr},  0  <  t  <  T. 

i£R  Jo 

Combining  (3.20),  (3.23)  and  (3.14),  we  find 

f  fj(k\w(x,t),v(x,t))dx  <  f  J7^(u>o(x), v0(x))dx  +  Ck2\ '^(0) 

J  —  OO  v  —  OO 

(3.25)  +  Ck  J  il/(r){  J  fj<'k\w(x,T),v(x,  r))dxjdT,  0  <  t  <  T. 

Then  the  Gronwall  inequality  yields  the  estimate 

/  fj(‘k)(w(x,t),v(x,t))dx 

J  —  OO 

<  |  J  fj{k)(w0(x),v0(x))dx 

(3.26)  +  CJb2y(fc)(0)|  exp(ckj  A/(r)dr),  0  <  t  <  T. 

Our  next  objective  is  to  take  the  £  power  of  relation  (3.26)  and  pass  to  the  limit 
k  — ►  oo,  in  order  to  obtain  a  relation  involving  the  L°°- norms  of  w  and  v.  To  this  end,  we 
first  estimate  the  right-hand  side  of  (3.26).  If  A m\z(w,v)  is  the  largest  eigenvalue  of  the 
Hessian  of  the  convex  function  f)(k)(w,v),  then 


(3.27) 


<  ^(  max  ^lx{w,v)){w2  +v2) 

l  V73  +  t;3<l 

<  \[  max  +  r?[*)(u;,r))](u;2  +  v2),  w2+u2<l; 

Z  u;3-fv3<l 


whence,  using  (3.15),  (3.16)  and  (3.19),  we  deduce 

(3.28)  fj{k)(w,v)<Ck2eCk(w2 +  v2),  w2  +  v2  <  1. 

Set  O  {x  £  R  :  “>o(x)  +  vo(x)  >  1},  and  observe  that  O  has  finite  Lebesgue  measure, 
0  <  [0\  <  oo.  Then,  using  (3.28),  (3.21),  (3.18),  (3.15)  and  (3.19)  we  arrive  at 

f  fj^k\w0{x),v0(x))dx  =  f  fj(-k)(w0(x),v0(x))dx 
J  —  oo  J  R\0 

+  I  fj{k\w0(x),v0(x))dx 

Jo 

<Ck2eck  f  (io2(x)  +  v2(x))dx 
J  -  OO 


(3.29) 


+  C|0|(1  +fc)K(fc)(0)sup  l;ie0(x)| 

ifR 

/wo(z)  _ 

vV(Orf£  '  }• 


Fi nally,  combining  (3.26)  with  (3.29)  and  (3.19), 


(3.30)  [  rj(k\w(x,t),v(x,t))dx  <  C(k2  +  k  +  1)  exp{k(A  +  (  M(t)cLt)}, 

J  —  oc  Jo 

where  the  constant  .4  only  depends  on  the  L°°-norm  of  the  initial  data  and  the  values 
of  if, '(w )  and  on  the  smallest  interval  where  ■u.’o(x)  and  uq  take  values.  Next,  we 

estimate  the  left-hand  side  of  (3.30).  We  set 


(3.31) 


and  consider  the  set 


5(f)  :=  sup { 1  +  | il' ( z ,  f ) |  +  ]v(z,f)|} 


E6(t)  :=  {x  £  R:  S(t)  -  S  <  1  +  |™(x,f)|  +  |v(x,f)|  <  5(f)}, 

where  S  >  0  small.  Then  | Es(t)j  <  oo.  Moreover,  by  virtue  of  (3.21),  (3.18),  (3.15)  and 
(3.17)  we  obtain 

f  Tj(k)(w(x,t),v(x,t))dx  >  j  Tj{k)(w(x,t),v(x,t))dx 
J— oo  J  E((t) 

>  f  T)(k\w(x ,t),v(x,t))dx 

Jem 

-  C(1  +  k)Yw{0)  f  (1  +  *)[)</* 

Jem 


(3.32) 


>  _  (7(i  +  k)Y{k){0)\E6{t)\S{t). 

C/ 


We  combine  (3.29)  with  (3.30)  and  raise  the  resulting  inequality  to  the  l/k  power. 
Then,  we  first  let  k  — *  oo  and  then  6  j,  0.  After  taking  the  logarithm  of  both  sides  we 
conclude 


(3.33) 


ip{l  -j-  |u;(x,f)|  +  |r(x,f)!}  <  C  -f  C  f  A/(r)dr,  0  <  f  <  T. 
--R  Jo 


Substituting  (3.24)  into  (3.33)  and  using  (3.31)  yields 


(3.34) 


<  C  +  C(1  +  T)  f  S(r)di 
Jo 


Finally,  integrating  (3.34)  we  conclude 


(3.35) 


sup{l  4-  |«7(x,  f)|  +  |r(x,  f)|}  <  A/,  0  <  t  <  T, 


where  the  constant  M  depends  on  T  but  not  on  c.  Therefore,  the  solution  {w(x,t),v(T,t)} 
of  (2.4)  is  bounded  on  ^0,  T]  x  R,  uniformly  in  e  >  0. 

The  a  priori  estimates  (3.14)  and  (3.35)  together  with  an  existence  theorem  for  the 
regularised  problem  (1.8)  are  summarized  below. 

Theorem  3.1.  Under  the  hypotheses  (1.5)-(l.7),  for  each  e  >  0,  T  >  0,  the  initial 
value  problem  (2-4)  (respectively  (1.8))  has  a  unique  solution  {ie(x, /),  r(x,  <)}  defined  on 
R  x  [0,7\  such  that  w,  v  £  C(t0,  Tj;  Z2(R))  n  Z,°°(R  x  [0,T]),  wx.vx,wxx,vxx,wt,vt  6 
C((0,  T};  Z2(R))  and.  also,  wx,vx  £  Z2(R  x  [0,Tj).  In  addition ,  as  e  \  0,  the  families 
{ic(x,  t),  r(x,  <)}£>u  and  {s^ir^x,  t),  s1/2vx(x,  t)}e>0  lie  in  bounded  sets  of 
X°°([0,T;;L2(R)  n  I*(R))  and  L2(R  x  [0,T]),  respectively. 

PROOF.  It  remains  to  prove  the  existence  part  of  Theorem  2.1. 

In  view  of  (3.35)  we  modify  the  function  as  follows.  Let  M  be  the  a-priori  Z—bound 
in  (3.35)  and  choose  i  >  max{M,  We  modify  the  function  <p( w )  outside  the  interval 

so  that  the  resulting  function  fi(w)  is  twice  continuously  differentiable  and  satisfies 


(3.36) 

and 

(3.37) 


-{ 


( p(w)  if  |tn|  <  £ 
a  linear  function  if  liwl  >  2 £ 


<  <p'(w)  <  1  +  max  ip'(w)  =:  L,  w  6  R. 

—  Kw<l 


Then  the  derivation  of  (3.35)  suggests  that  any  solution  {u;(x,  <),  u(x, /)}  of  the  system 
(2.4)  with  ip  replaced  by  will  satisfy  the  bound  (3.35)  with  the  same  constant  M.  (We 
note  here  that,  although  the  family  of  entropies  for  the  modified  system  (2.4)  differs  from 
the  original  (cf.  (3.15),  (3.16)),  the  estimates  (3.17),  (3.18),  (3.19),  (3.29)  and  (3.32)  hold 
with  the  same  constants  for  both  families).  By  virtue  of  the  above  remarks  we  may,  without 
loss  of  generality,  assume  that  the  function  <p(w)  in  (2.4)  has  been  modified  to  comply  with 
(3.36)  and  (3.37). 

Existence  and  uniqueness  of  solutions  for  (2.4)  on  Qt  =  Rx  [0,  T]  will  be  deduced  by 
a  simple  application  of  the  contraction  mapping  theorem.  Let  B  be  the  Banach  space 

B  :=  {(u.-,r)  :  Or-  -  R  x  R,u7,r  6  Z°°([0, T*]; Z2(R)), 
wt,vx  6  Z2([0,T*];Z2(R))} 
for  some  0  <  T*  <  T,  with  norm 

/oo  »X*  <>oo 

[w2(x,t)  +  v2(x,t)}dx  +  e  /  /  (w2x  +  v\)dxdt. 

_  _  -oo  */ 0  */  —  oo 

Let  S  be  the  map  that  carries  (IF,  V’)  £  B  into  the  solutions  (w,v)  of  the  linear  parabolic 
initial  value  problem: 

’  wt  =  VX  +  ewxx 

(3.38)  <  vt  =  <f( W)x  +  T\V}  + evxx 

u’(x,0)  =  U’o(x),  v(x,0)  =  t’o(x). 


13 


By  standard  theory  of  the  heat  equation  S  :  B  — ♦  B,  for  any  T*  >  0.  We  proceed  to  show 
that,  for  T*  small  enough,  5  is  a  contraction. 

Let  c-  b  and  let  (wj ,  vj ),  (w2 ,  U2)  be  the  respective  solutions  of 

(3.38)  corresponding  to  the  same  initial  data  itq(2,0)  =  u>0(i),  Vj(2,0)  =  ^0(2),  i  =  1,2. 
Simple  energy  estimates  for  (3.38)  yield 

/OO  ft  fOO 

(u’1(2,f) -W2{x,t))2dx  +  £  /  /  (wi(x,t)  -  w2(x,T))ldxdT 

-  00  J  0  J —  00 

(3.39)  <  -  T  [°°  (\\(x,r)-V2(x,T))2dxdT, 

£  J  0  J  —  OO 


/00  W  *oo 

(t>l(X>0  -  1-2(2,  t))2dx  +  £  /  /  (Vi(x,T)  -V2{x,T))ldxdT 

-  OO  J 0  V  —00 

<  ~  f  I"  W{W\{x,T))-tp{\V2{x,T)))2dxdT 
€  J  0  J  —  00 

+  {v!(x,t)  ~v2(x,T))2dxdr 

J0  J  —  00 

+  f  r  (r[V2(x,r)-V2(x,r)])2dxdr. 

JO  J  —  00 


(3.40) 


Moreover,  in  virtue  of  (3.3 7) 


(3.41) 


and  by  (1.4) 


r°° 

/  /  (^(^(2,-r))  -  v?(W2(2,r)))2d2dr 

«/o  J —00 

<  L2  f  f  (\\\(x,t)  -  W2(x,t))2 dxdr 
Jo  J- 00 


/  /  (^[VHx.t)  -  V2(2,r)])2d2dT 

(3.42)  <2[max(|r(t)|  +  |r'(<)|)2](r2  +  l)  f*  H  (Vy(x,r)  -  V2(x,t))2 dxdr. 

[o,T]  y0  «/ — 00 

Combining  (3.39)-(3.42)  and  using  the  definition  of  norm  in  B ,  one  obtains 

lll((»i  -  «-,),(»,  -  >^)):i!i(0  <  c.r-i!!((ir,  -  ir2),(v.  -  v,))!!&(n 

(3.43)  +  f  |||((t»,  -tlli),  (v,  - 

Jo 


i»!V?S«!vr 


By  the  Gronwall 


£ 


> 

>> 
> 


;■  a 

y.v 


v  - 

Cv!- 


V’ 


£ 


.V. 


>'a‘ 

.  A' 

■v. 


rr. 


where  C£  >  0  is  a  constant  depending  on  ^ ,  X,  X,  and  !|-ft||c1[o,T’J- 
inequality, 


(3-44) 


III ((™1  -  -  W2))liifl(T“) 

<ceT*i|i((iri-u'2),(r1-v2))i!!2B(T*), 


and  S  is  a  strict  contraction  on  B,  provided  X*  is  chosen  sufficiently  small. 

By  proceeding  in  steps  of  length  X*,  the  above  procedure  yields  a  unique  solution 
{u;(x,t),v(x,  f)}  on  R  x  [0,  Xj  of  the  system  (2.4),  for  each  fixed  e  >  0  and  any  X  >  0.  By 
standard  X2-theor y  of  the  heat  equation  [10],  the  solution  w(x,  t),v(x,  t )  £  C(j0,  X];  X2(R)) 
and  has  the  regularity 


(3.45) 

(3.46) 


Vi  wx,  V~tvx£  X°°([0,Xj;  X2(R)), 
wx,vz,Vt  wt,Vt  vt,Vt  wxx,Vt  vxx  £  X2(R  x  [0,  X]). 


The  solution  {ie,u}  enjoys  additional  regularity  which  is  obtained  as  follows.  Differ¬ 
entiate  (2.4)i  and  (2.4)2  with  respect  to  t  (the  rigorous  justification  is  by  taking  difference 
quotients)  and  use  (1.4)  to  obtain 


* 

i 

(3.47) 

35!  • 

£ 

(3.48) 

wtt  =  vxt  +  ewxxt, 


vtt  =  <e(u')xt  4-  r(0)ut(x,  t)  +  /  r’(t  -  r)u((x,r)dr  +  evxxt. 

Jo 


Multiply  (3.47)  by  t2wt  and  integrate  over  R  X  [0,X|.  Integration  by  parts,  (3.35),  (3.45), 
(3.46)  yield 


/OC  fOO 

u?2(x,  t)dx  +  /  /  T2w\tdxdr 

-oo  JO  J  —  oo 

(3.49)  <  Ce  f  f  ( rv 2  +  TW2)dxd,T  <  Cc. 

Jo  J -oo 

Similarly,  multiplying  (3.48)  by  t2vt  and  performing  the  same  steps,  we  obtain 

/oo  ft  too 

v2(x,t)dx  +  I  I  r2v\tdxdr 

-oo  J  0  J  —  oo 

(3.50)  <  Ce  (  f  ( tv 2  +  rw\)dxdr  <  Ct. 

Jo  J- oo 
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Next,  we  multiply  (3.47)  and  (3.48)  by  i3wtt  and  t3 rtt,  respectively,  integrate  by  parts  over 
R  x  [0,  t]  and  use  (3.35),  (3.45),  (3.46),  (3.49),  (3.50)  and  the  Cauchy-Schwarz  inequality 
to  arrive  at 

A  t  A  OO  m  OO 

I  /  T3w\tdxdr-\-t3  j  wlt(x,t)dx 

J  0  J  - oo  J  —  oo 

(3.51)  <Ce  f  f  {r2v2xt  +  T2w\t)dxdT  <  Ce, 

JO  J  -  oo 


J 


(3.52) 


Al  AOC  aQO 

I  j  T3v2tdxdr  +  t3  /  v2xt(x,t)dx 
Jo  J  —  oo  J  —  oo 

<  Ce  f  f  (r2v2t  +  r3(iyp(w)xt)2  +  tv2 )dxdr 

Jo  J  —  OO 

<  Ce(  1  +  f  (  T3wlw2dxdr) 

JO  J~co 

<  Ce(l  +  f  rsup  w2x(x,  r)dr) 

Jo  * 

<  Cc(  1  +  f  t  f  (wl  +  wxx)dxdr)  <  Cc. 

JO  J  —  oo 


Finally,  using  (3.45),  (3.46),  (3.49),  (3.50),  (3.51),  (3.52)  and  standard  embedding  theorems 
we  conclude  that  {u;(x,  t),  v(x,  <)}  has  the  regularity 

claimed  in  Theorem  3.1. 

4.  THE  PASSAGE  TO  THE  LIMIT  e\  0. 

For  s  >  0,  let  {we(x,  t),  vc(x,  t)}  be  the  solution  of  the  initial  value  problem  (1.8)  on 
Qt  :=  R  x  [0,T],  with  regularity  properties  as  in  Theorem  3.1.  By  virtue  of  (3.35)  the 
family  of  functions  {wc(x,  t),ve(x,t)}£>o  is  uniformly  bounded  in  L°°(Qt)-  There  exist 
functions  w(x,t),v(x,t)  and  tp(x,t)  in  L°°(Qt)  such  that,  along  a  subsequence, 

(4.1)  we(x,t )  w(x,t )  in  L°°  —  weak  star, 

(4.2)  ve(x,i)  u(x,t)  in  L°°  -  weak  star, 


(4.3)  tp(we(x,t))  <f(x,t)  in  L°°  -  weak  star, 

as  e  \  0.  The  object  of  this  section  is  to  show  that  {u;(x,  <),  r(x,  <)}  is  a  solution  of  (1.1) 
in  the  sense  of  distributions. 


iW 


The  major  obstacle  to  overcome  is  that,  in  general,  nonlinear  functions  are  not  con¬ 
tinuous  under  weak  star  convergence,  and  it  does  not  follow  that  i p(x,t)  =  <p{w(x,t)). 
Under  weak-star  convergence,  such  composite  weak  limits  have  been  characterized  bv  Tar¬ 
tar  [21,22,23’  as  expected  values  of  a  family  of  probability  measures  zqIt),  called  the 
Young  measures.  The  relevant  issue  here  is  whether  the  Young  measure  reduces  to  a  Dirac 
mass.  Using  the  theory  of  compensated  compactness  and  a  class  of  entropy-entropy  flux 
pairs  introduced  by  Lax  ill],  DiPerna  [7]  proved  the  following  key  resuit  which  will  be 
employed  in  the  sequel. 

Proposition  4.1  (DiPerna).  Let  {u>e,uE}  :  Qj  — >  R  be  a  collection  of  functions  such  that 

!i’u,e  !Iloo(<?t)  +  !ke|!z.“(QT)  ^ 

where  C  is  a  constant  independent  of  e.  Suppose  also  that  for  any  smooth  entropy- entropy 
flux  pair  {r](w,  v),  q(w,  u)}  of  (E)  (see  Section  3),  with  <p  satisfying  (1.5), 

dtTj(we(x,t),ve(x,t))  +  dxq(w‘(x,t),ve(x,t )) 

lies  in  a  compact  set  of  H  i0]{Qt)-  Then  there  exists  a  subsequence  {w~  ,vc  }  and  functions 
w,v  £  L°°(Qx)  such  that 

wc  ( x,t )  — >  w(x,t),  v£  (x,t)  — *•  v(x,t),  a.e.  for  (x,  t)  S  Qr 

as  e'  l  0. 

We  apply  Proposition  4.1  to  the  family  {we(i,  t),ve(x,  <)}e>o  of  solutions  of  (1.8) 
(which  also  satisfy  (2.4)).  A  straightforward  calculation  using  (2.4)  and  (3.2)  yields 

dtT)(wc,vc)  +  dxq{wc,vc ) 

=  e^2dx(ei/2Vw(wty)wez  +  e1/2r?vK,ueK) 

-  £[f?u>w(ttfe,t’*)(w*)2  +  2 rjwv(we  ,ve)wexvex  +  T)vv(we ,v')(vex)2} 

(4.4)  4-  T)v(we  ,ve)T\vc). 

The  a-priori  estimates  (3.14)  and  (3.35)  imply  that  the  family 

51/23«(e1'2ii 1v,{w',v')wi  +  e^27fv(wc,ve)v‘) 

converges  to  0  and  is  thereby  compact  in  while  the  family 

e{t]ww(we  ,vc)(w£x)2  +  2riwy(u)€  ,v')wcxvx  +  77vv(u,’e,u£)(v£)2] 
is  bounded  in  Li(Qx)-  In  addition,  the  family 

(4.5)  dtj i(we,ve)  -f  dzq(we,ve)  -  qv(we,  u£).F[Vj 


resides  in  a  bounded  set  of  H'~1,00(Qt)-  Using  a  lemma  of  Murat  [14]  we  deduce  that  the 
quantity  (4.5)  lies  in  a  compact  set  of  Finally,  since  rjv(we  ,vc)T{v‘\  lies  in  a 

bounded  set  of  L2{Qj)  and  thus  in  a  compact  set  of  H[q1c(Qt),  the  family 

dtT)(wc,ve)  +  dxq{wc,vc) 

also  lies  in  a  compact  set  of  HfJffQr)-  Thus,  Proposition  4.1  implies  that 

we  ( x,t )  — >  w(x,t),  ve  (z,t)  — >  v(x,t)  a.e.  in  Qt, 

along  a  subsequence  e'  \  0,  and  permits  passage  to  the  limit  s'  \  0  in  (1.8)  in  the  sense 
of  distributions.  The  pair  of  functions  {w(x,  t),  v(x,  t)}  belongs  to  L°°([0,  T];  L2(R))  n 
L°°{Qr)  and  is  a  weak  solution  of  (1.1)  on  Qt-  This  completes  the  proof  of  Theorem  1.1. 

Finally,  we  remark  that  one  does  not  expect  weak  solutions  of  (1.3)  to  be  unique.  We 
note  that  for  every  (convex)  entropy  r)(w,v),  the  solution  constructed  here  satisfies  the 
inequality 

dtT)(w,v )  +  dxq(w,v)  <  qv(w,v)T{v\ 
in  the  sense  of  distributions. 
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